Theoretical analysis of the transmission phase shift of a quantum dot in the presence 

of Kondo correlations 
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We study the effects of Kondo correfations on the transmission phase shift of a quantum dot 
coupled to two leads in comparison with the experimental determinations made by Aharonov-Bohm 
(AB) quantum interferometry. We propose here a theoretical interpretation of these results based 
on scattering theory combined with Bethe ansatz calculations. We show that there is a factor of 2 
difference between the phase of the S-matrix responsible for the shift in the AB oscillations, and the 
one controlling the conductance. Quantitative agreement is obtained with experimental results for 
two different values of the coupling to the leads. 
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Quantum dots (QD), small puddles of electrons con- 
nected to leads, can be obtained in a controlled fashion 
thanks to recent progress in nanolithography. Under cer- 
tain conditions a dot can be modeled as a localized spin 
coupled to Fermi baths (the leads). A Kondo effect takes 
placoii^^ when the temperature is lowered . A key ingre- 
dient of the Kondo effect is the phase shift 8 an electron 
undergoes when it crosses the dot. While its direct mea- 
surement was out of scope in bulk systems, it became fea- 
sible recently in quantum dots via Aharonov-Bohm (AB) 
intcrferometery4. We mention here the experimental re- 
sults obtained in two cases corresponding to a strong cou- 
pling to the leads^. In the unitary limit, the phase shift 
climbs almost linearly with Vq with a value at the middle 
of the Kondo enhanced valley which is almost ir. At a 
smaller value of the coupling strength, the phase shift de- 
velops a wide plateau at almost tt. We will call the latter 
case the " Kondo regime" . Early theoretical work on the 
phase shift for the bulk Kondo effeci&£ predicts 6 = tt/2. 
In the context of QD, Gerland et a& had obtained, on 
the basis of NRG and Bethe ansatz calculations, a vari- 
ation of 8 with the energy of the localized state leading 
to a value of tt/2 in the symmetric limit, in disagreement 
with the recent experimental results quoted above^^. In 
this paper, we propose a new theoretical interpretation of 
the experimental results based on scattering theory and 
Bethe ansatz calculations. Our main prediction concerns 
a factor of 2 difference found between the phase of the 
S-matrix observed by the phase shift measurements and 
the phase governing the conductance. 

Let us consider a quantum dot coupled via tunnel 
barriers to two leads L and R, and capacitively to a 
gate maintained at the voltage Vq. The system can 
be describedSiifi by an one-dimensional Anderson model 
with two reservoirs L and R 
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Consider the elastic component of the S-matrix, Ska-, 
describing the scattering of a spin-tr electron with mo- 
mentum k off the impurity. It is given bj»iii£iifi Sk 

CM 



>fc<T — 

iT^ s ), where C CT is a multiplicative phase factor 
and T^ s is the T-matrix with matrix elements given by 



(2) 



where a, (3 — L or R, p a (sk) is the density of states of 
conduction electrons for a and and G a ( £ k + i-v) is 
the exact localized electron retarded Green's function. 
Using exact results&iA^ on the self-energy at T = in 
an interacting Fermi liquid, one can show that riQ a = 
^Im\nQ a (fj, + irj). Friedcl's sum rulei 3 .^ requires no a 
to be equal to the change in the number of conduction 
electrons with spin a resulting from the addition of the 
impurity. Hence it is related to the transmission phase 
shift 8 a at the Fermi level, no CT —— ~ 8& . Therefore 8 a coin- 
cides with the phase of the Green's function at the Fermi 
level Q a (n+irj) . If we denote the associated self-energy by 
T llJ (fi+ir]), one gets Q a (n+irj) = sin S^e 15 " /ImT, cr (fi+irj). 
with 7mS CT (/i + ir,) = -tt(U l 2 + V%)p a {^^ at T = 
leading to 
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In the case of a symmetric QD with Vt, = Vr 



has 



Cni sin 8„ 



and S^ a 
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= V, one 
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C CT cos 8 c % ^ a . The multiplicative phase factor C a con- 
tains additional information about the spectrum and the 
filling of the quantum dot. To determine it, we will make 
use of Levinson's theoremi4iA£. In its original form, the 
theorem applies to the potential scattering of a parti- 
cle in a given momentum I and relates the zero-energy 
phase shift 8i to the number of bound states of the same 
I supported by the potential. It was gcncralizedi&iii later 
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on to the case of the scattering of a particle by a neutral 
compound system as constituted for instance by an atom. 
In the present case of a QD which can be viewed as an 
artificial atom, it follows that lndet S , fc J7IT /(2i7r) is equal 
to the total number of states i.e. no a = uq. By ap- 
plying generalized Levinson's theorem to Sk F a, one finds 



a 



and 



S, 



cos Scr i sin 5 a 
i sin Srr cos 5„ 



(4) 



where 5 = 5 a . One can easily check that, Sk FIJ be- 
ing a unitary matrix, the optical theorem is fulfilled: 
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FIG. 1: Experimental conductance G e x P (Va) and phase shift 
v(Vg) as a function of Vg (values taken from Ref0 cf. text). 
Comparison is made with the curve G(Vg) = sin 2 (^j(Vg) / '2). 

In an open Aharonov-Bohm interferometry 
experiment^, spin-tr electrons coming from the source 
through each of the two arms interfere coherently 
at the drain, leading to periodic oscillations of the 
differential conductance, the argument of which is 
given by 2n$e/h + 5qo- $ is the magnetic flux and 
Sqd is the transmission phase shift introduced by the 
QD, equal to 5 = ttuq (cf. Eq0}. In this paper, we 
neglect the role of the reference arm on the phase shift 
considered by some authors^ and concentrate on the 
contribution of the quantum dot to the interference 
pattern. The conductance through the QD is expressed 
by the Landaucr formulaic, G oc J2a \ T k*J 2 - Usin g 
EqQ] we get G oc sin 2 5 a . In the absence of magnetic 



field, <5| = <5j = 5/2, one gets 

Goc ^sin 2 5/2. 



(5) 



Due to recent developments in experimental tech- 
niques, one now disposes of simultaneous measurements 
of G and 5qd- In this paper we check the validity of 
the theoretical prediction of Eas. l|4l5|l by reporting the 
experimental results for G and 5qd obtained in the uni- 
tary limit at different values of Vg- Before examining the 
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FIG. 2: (a) Theoretical results from the B.A. calculations for 
the occupation number no as a function of the normalized 
energy (so/U + 1/2) at different values of T/U. Note that 
no = 1 at the symmetric limit eo = —U/2 and the existence 
of a plateau in its vicinity when T/U < 0.25; (b) The same 
quantity as a function of the renormalized energy e* /Y at 
different values of F/U. 



experimental test, we make the following remarks: (i) in 
an interferometry experiment, only relative values of the 
transmission phase shifts can be measured. Hence we set 
5 = 7T at the location of the maximum of the visibility, 
evaluated to Vg = 423mV. This implies a shift in the 
(5-scale evaluated to A5 = 0.297T with ip = 5 + A5; (ii) 
typically the measurement of the conductance G is per- 
formed in a "one-arm" device (pinching off the reference 
arm with the barrier gate) whereas that of the visibil- 
ity is done in a "two-arm" device. As a result, while 
the evolution of the visibility with Vg mimics that of the 
conductance, the value of the former is shifted with re- 
spect to that of the latter, by AVg = 15mV. Therefore 
we take the values of G at (Vg — AVg), and of 5 at Vg; 
(iii) the conductance is normalized by its maximum value 
at Vg = 423mV. Taking all these points into account, 
the graph reported in Fig^ shows that the experimental 
dependence of sin 2 <p/2 with Vg reproduces that of the 
"shifted" conductance G exp in a quite remarkable way, 
providing further support to the validity of Eqs.003^- 
We now want to evaluate no in order to derive 5 = imo . 
Starting from Eq. Q , one can show^ that only the sym- 
metric linear combination of electrons couples to the lo- 
calized state. Therefore if we are only interested in no, 
it is sufficient to study a single reservoir Anderson model 

We 



with a hybridization potential V = y/V^ + V^ 
have solved the equations of the Bethe ansatz (B.A.) 
numerically at T = 22 i 23 i 24 . This allows us to deter- 
mine the value of no as a function of the parameters of 
the Anderson model e , V and U . The three param- 
eters enter through their ratios Eq/U and T/U, where 
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FIG. 3: Fit of the experimental data for the VG-dependence of 
the phase shift with the B.A. results for 5 = nno as a function 
of — Eg/r. Making use of the electron-hole symmetry, experi- 
mental results both below and above the symmetric limit are 
reported using the same scales. They are represented by tri- 
angles pointing down and up respectively in the unitary limit, 
and by circles and squares respectively in the Kondo regime 
(incorporating a shift in the 5-scale cf.text). The best fit is 
obtained for T/U = 0.5 in the unitary limit (both below and 
above the symmetric limit), and for T/U = 0.07 or 0.05 in 
the Kondo regime (below or above the symmetric limit re- 
spectively) . 
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FIG. 4: Phase shift as a function of the gate voltage Vg- 

(a) Unitary limit. Theoretical results from Bethe ansatz cal- 
culations at T/U = 0.5 compared to the experimental data 
for tp/n = 8 exp /ir + 0.29 (triangles pointing down and up) 

(b) Kondo regime. Same thing with T/U = 0.04 and 0.07, 
Lp/n — S exp /-K + 0.01 (circles and squares). 



r = ttV 2 pq. Denote by no(eo,T /U) the value of no for 
the corresponding values of the parameters. The follow- 
ing relation holds due to the particle-hole symmetry of 
the model^: n (-(e + U),T /U) = 2- n Q {e Q ,T/U). 
This automatically ensures no(—U/2,T/U) = 1 in the 
symmetric limit So = —U/2 whatever T/U is. Further- 
more it follows from the preceding relation that the study 
can be restricted to —U/2 < Eq < U/2 and the remain- 
ing part can be deduced from it. The results of the 
calculations are reported in Fig. 2(a). For strong cou- 
pling strengths T/U > 0.25, no is found to climb almost 
linearly with ~(eq/U + 1/2) whereas for weak coupling 
strengths T/U < 0.25, the energy dependence of no de- 
velops a plateau around Eo = —U/2. This change of be- 
havior is due to the fact that the extent of the local mo- 
ment regime (centered around eo = —U/2 with no — 1) 
increases when T/U decreases. As the temperature is 
lowered, the Kondo resonance develops through this lo- 
cal moment regime. This plateau-like structure can be 
viewed as the beginnings of the "staircase" variation of 
no with (To obtained in the localized regime T/U — * 0. 

The experimental data can be fitted then with two pa- 
rameters, T/U and Eq/U. The value of eo is governed by 
the strength of the gate voltage. Fitting the experimen- 
tal data from results presented in Fig. 2(a) is a difficult 
task since one needs to fix the correspondence between 
£o/U and Vg on the one hand (we take it linear as usual, 
independent of the regime considered), and to find the 
best fitting value for T/U in the different regimes on the 
other hand. A valuable help for doing this is provided 
by taking advantage of some special properties of the 
Anderson model. These properties can be easily recog- 
nized when physical quantities such as no arc plotted 
as a function of some renormalizcd energy defined as 
e* /T = £q/T + g(U/(T)). In the asymmetric regime^ 



when (U + 2e ) > vTTJ and |e | < U, g{U/{T) equals 
^ ln(7re£//(4r)) and the behavior of no as a function of 
— £q/T is universal^^i. This property is illustrated in 
Fig. 2(b). The universality is reached when T/U < 0.25 
and the range of energy over which universal behavior 
extends is given by |eg/r - 1/tt ln(aU/T) | < U/T. One 
can also see from Fig. 2(b) that in the empty level regime 
(?io — ► 0), the curves n = /(eq/T, T/U) at various values 
of T/U merge, displaying an asymptotic behavior—. The 
existence of both these universal and asymptotic behav- 
iors is of valuable help in fitting the experimental data. 
FigCHrcports the results of the fit in the unitary limit and 
Kondo regimes. The experimental results incorporate a 
shift in the 5-scale, ip — S + A5 in order to get cp = tt 
at the symmetric limit. We establish the correspondence 
between Vg and £q/U by fitting the experimental data to 
the theoretical results in the empty level regime when all 
the curves merge. One finds AVg/A(£o/l^) of the order 
of 30mV in both of the regimes considered. The best fit 
is obtained for T/U = 0.5 in the unitary limit both below 
and above the symmetric limit, and for T/U = 0.04 or 
0.07 in the Kondo regime (below or above the symmetric 
limit respectively). Finally by keeping the same corre- 
spondence between Vg and £q/U and using <5 = 7rno, we 
derive the dependence of the phase shift with Vg from re- 
sults obtained in Fig. 2(a). As can been seen from Fig0J 
our theoretical predictions are in quantitative agreement 
with the experimental data. The fit is all the more re- 
markable that it is performed in presence of a single fit- 
ting parameter T/U only. 

In conclusion, we have proposed a theoretical analy- 
sis of the transmission phase shift of a quantum dot in 
presence of Kondo correlations and confronted our re- 
sults with the Aharonov-Bohm intcrfcrometry and con- 
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ductance measurements. We have explained the presence 
of a factor of 2 difference between the total phase of the S- 
matrix (responsible for the shift in the A-B oscillations), 
and the one appearing in the expression of the conduc- 
tance G ~ sin 2 (<5/2). Our calculations based on Bethe 
ansatz lead to a remarkable quantitative agreement with 
experimental results. The whole discussion so far has 
been restricted to the low temperature regime. One of 
the next goals will be to include finite temperature ef- 
fects as well as to study the role of a magnetic field and 
consider the out-of-equilibrium situation. 
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